Heat conductivity from molecular chaos hypothesis in locally confined billiard systems 
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We study the transport properties of a large class of locally confined Hamiltonian systems, in which 
neighboring particles interact through hard core elastic collisions. When these collisions become 
rare and the systems large, we derive a Boltzmann-like equation for the evolution of the probability 
densities. We solve this equation in the linear regime and compute the heat conductivity from a 
Green-Kubo formula. The validity of our approach is demonstated by comparing our predictions 
to the results of numerical simulations performed on a new class of high-dimensional defocusing 
chaotic billiards. 
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The understanding of Fourier's law and the computa- 
tion of the heat conductivity in Hamiltonian systems as 
a function of temperature and of the physical parameters 
remains to this day a challenging issue [lj. In particular, 
establishing the necessary conditions the dynamics must 
satisfy so as to justify a first principles based derivation of 
Fourier's law has been the subject of ongoing discussions. 

As a generic model of heat transfer in insulating crys- 
talline solids, one often considers a lattice of coupled par- 
ticles with nearest-neighbor interactions whose motion 
obeys Hamilton's equations. Thus consider TV particles 
of unit masses located on a one-dimensional lattice with 
positions and momenta (q, p) = {(qi, Pi)} 1<i <j V ' w ^ n 
q;, pi 6 R d . The Hamiltonian H takes the form 
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where V represents the interaction with the external sub- 
strate and U the nearest- neighbor interactions [Toj ] . 

After Peierls' work 0, all attempts to give a satisfac- 
tory derivation of Fourier's law in mechanical systems 
have focused on the study of weakly anharmonic dynam- 
ics. Using the Peierls-Boltzmann equation, recent works 
have studied the effects of phonon collisions on the heat 
conductivity |, |, i, @, . In this context, the conductiv- 
ity may be interpreted as a collision frequency between 
phonons. 

In this letter, we focus on the opposite limit, namely 
extremely anharmonic interactions, and, under minimal 
assumptions on the chaotic nature of the dynamics, iden- 
tify a class of models which display a universal response 
to non-equilibrium thermal constraints. The motivation 
for this study is twofold: First, the heat conductivity 
can be computed from first principles and takes a sim- 
ple form; Second, as pointed out in |8j, such systems of 
locally confined particles in interaction find concrete ap- 
plications in the study of aerogels, materials in which 



gas particles are trapped in nano-size pores and rarely 
interact among themselves. Assuming the validity of 
a Boltzmann-like equation to describe such systems of 
rarely interacting particles when they become large, we 
show that the heat conductivity of such systems is gener- 
ically equal to the average frequency of interaction be- 
tween the systems' components, i.e. irrespective of the 
detailed geometric properties of the confinement mecha- 
nism. This will be checked in detail by numerical simu- 
lations, showing the universality and power of the Boltz- 
mann approach to analyze the transfer of heat in the 
mechanical systems we study. 

To be specific, we consider the case of interaction po- 
tentials which take only the values zero inside a region 
Vi,u C M. d with smooth boundary A of dimension d — 1, 
and infinity outside. Likewise, the pinning potential V is 
assumed to be zero inside a bounded region fly and infin- 
ity outside, implying that the motion of a single particle 
remains confined for all times. The regions flu and Qy 
being specified, the dynamics is equivalent to a billiard in 
higher dimension. An important quantity in such models 
is the average rate of collisions between nearest-neighbors 
under equilibrium conditions. We will be specifically con- 
cerned with the limit of rare collision events. 

The shape of the region f2v_ determines the nature of 
the local dynamics. In ref. 8], fly was chosen to be a 
semi-dispersive billiard with bounded horizon, thus en- 
suring strong chaotic properties of the dynamics. In par- 
ticular the fast decay of correlations of the local dynam- 
ics was invoked to set up a stochastic equation describing 
the energy exchange dynamics. It is our purpose to show 
that this assumption can be relaxed: local ergodicity is 
enough to warrant the identity between heat conductivity 
and frequency of energy exchanges. We regard this as an 
important result which further validates the analogy be- 
tween this class of models and aerogels whose nanopores 
need not have dispersing properties. 

Examples of the simplest type of billiards we may con- 
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sider are periodic arrays of square boxes in two dimen- 
sions in each of which a single hard disk particle moves 
freely, but can still perform collisions with neighboring 
disks by interacting through the confining walls, for in- 
stance, provided we let the cells overlap a bit. The 
specific nature of the interaction mechanism at play is 
however not relevant in our formalism. We will instead 
consider point particles moving freely in two-dimensional 
square boxes of unit sides and interacting among near- 
est neighbors when the Euclidean distance between them 
becomes equal to a parameter which we denote by a. At 
that point, they exchange their longitudinal velocities, 
i.e. the velocity components in the direction of their rela- 
tive motion. We refer to this model as the square-strings 
model. The interaction may be depicted by attaching 
strings of lengths a separating neighboring particles, as 
shown in Fig. [I] In this case, we take Qy — [—1/2, 1/2] 2 , 
and flu = P? _i n)( a )i the disk of radius a with center 
at (—1,0) We note that, in the absence of interac- 
tions, the dynamics of the individual particles is pseudo- 
integrable; it is ergodic on the configuration space for 
most values of the velocity directions, but is known to be 
non-mixing. We will consider this model in some details 
below and provide numerical evidence that the analysis 
which follows applies to it. 
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FIG. 1: (Color online) Typical trajectories of the square- 
strings model, displayed for increasing values of a, color-coded 
from blue to red according to their energies. Small and large 
arrows indicate initial and final velocities respectively. For 
large a, the system appears to be near-integrable, reflecting 
the rarity of interactions, but is nevertheless fully chaotic. 

This model can be compared to a simpler class of com- 
plete exchange models, specified by square well poten- 
tials, obtained from Eq. ([1]) with d = 1, as a limit of 
models with smooth interaction potentials. In this case, 
Civ = [—b,b] and fljj — [—a, a]. Each particle on the 
lattice moves freely on a one-dimensional cell of size 2b, 
changing directions at the boundaries. The interaction 
between a pair of particles acts when the difference be- 
tween the positions of the two particles reaches the value 
a, at which point they exchange their velocities. 

In the general d-dimensional set-up, the particles move 



freely inside their respective cells, bouncing off the walls 
elastically, until the vector q; — q i+ i (resp. q^_i — q^) 
reaches the boundary A. The corresponding particles 
then exchange the components of their velocities in the 
direction normal to the boundary A, i.e. longitudinal to 
the direction of their relative motion. 

The Hamiltonian ([1]) may be written as a 
sum of local terms, hi(p, q) = pf/2 + V{c\j) + 
1/2 [U (qi_i - qi) + U (q s ; - q i+ i)] . This allows one 
to define a function describing the local transfer of 
energy by computing the variation in time of the 
local energy hi along the solutions of the equations 
of motion, 4;hi(p,q) — ji-i — ji, where the local 
energy current between sites i and i + 1 is defined as, 
ji = h(pi + Pi+i) ■ Vt/(qj - q 4 +i), which, for hard core 
interactions, becomes 

ji = -^ A (q< - q i+ i)brVi - [{Phi? - (pi) 2 ] , 

(2) 

where \x\ + = x, if x > 0, and otherwise, pf- = • n 
is the component of the vector p; in the direction of the 
unit vector n, normal to the boundary A, and S\ denotes 
the delta function concentrated on this boundary. The 
first factor corresponds to the localization of the collisions 
in configuration space, the second one gives the rate at 
which collisions occur and the last one corresponds to 
the exchange of longitudinal components of the kinetic 
energies. 

Starting from the Liouvillc equation for the evolution 
of probability densities on phase space, it is straightfor- 
ward to derive an equation for the evolution of the proba- 
bility density of a single particle in a given cell. It involves 
the probability distribution of the pairs of particles which 
consist of the particle itself and either of its nearest- 
neighbors on the lattice. The Boltzmann approximation 
simply amounts to assuming that this two-particle distri- 
bution factorizes in terms of the one-particle distributions 
fi at each site. To justify this assumption, one needs to 
show that a version of molecular chaos holds in our mod- 
els. Namely, that the dynamical variables involved in 
the successive collisions between two neighbors are in- 
dependent at the times of collisions. For that purpose, 
we require two ingredients: First, local correlations are 
typically destroyed after a collision between neighboring 
particles; Second, the number of particles must be very 
large, so that in the long run, the whole system plays the 
role of a reservoir for the specified pair of nearest neigh- 
bors. How these conditions are realized in the models we 
consider and, in particular, in the square-strings model 
which we test numerically, is not yet fully elucidated. We 
interpret the first condition as requiring interactions to 
be rare compared to the collisions within a single cell. In 
the square-strings model, it amounts to taking the maxi- 
mal separation close to the length of the diagonal joining 
opposite corners of neighboring boxes (= \/5), as in the 
third panel of Fig. [1] 
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We denote by / = {/j(p,q,i)} 1<i<JV , the set of the 
marginal probability distributions of each particle in each 
cell. The Boltzmann equation for this set of probability 
densities is 



dt 



fi(p, q, t) = -P • V q /, + V»U + Lt l+1 (/) + Lli_ x {f). 

(3) 

Here L w accounts for the collisions of the particles with 
the walls of their respective cells, and L^ i±1 for the col- 
lisions of the i-th. particle with the i ± 1th, viz. 

LU±i(L) = /d Pad q'S A (q- q')^ - p^\ + (4) 

x[/i±i(p&,q')/i(Pc>q) - /i(p,q)/i±i(p ,q')]) 

with p£ = p ± , p^ = p£, p c - pjrn = p p^n, and 
Pfc— p^ri = p a — pfn. One can check that the distribution 
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Mo q sIJ/i(Pi,q<) = ^- i n e_/3pr/2l ^(*) (5) 

1=1 1=1 

is stationary for any inverse temperature (3. Applied to 
this distribution, the advection term in Eq. ([3]) is zero 
except on the cell borders where it cancels with L w fi. 
(3 may be fixed by imposing identical thermal boundary 
conditions at both ends of the lattice. 

When the system is set out of equilibrium by imposing 
different temperatures at its boundaries, we proceed with 
a standard Chapman-Enskog expansion around a local 
equilibrium distribution, 

N N 

Wcq = J] /fc(p fc ,q*) = Z' 1 [] e-^/ 2 l Qv (q k ), (6) 

k=l fc=l 

with (3k — f3(k/N) for some smooth function (3, tak- 
ing as a small parameter the local temperature gradient. 
Plugging Eq. ([6]) into ([3]), we observe that only terms of 
second-order in the temperature gradient survive. This 
is in contrast with the case of an ordinary gas of colliding 
particles. This simplification occurs because the advec- 
tion term of the Boltzmann equation ([3]) acts only on the 
position variable within each cell and therefore not as a 
gradient on the lattice dependent variables. This means 
that local averages with respect to the distribution (J6j> 
are identical to local averages with respect to the true 
non-equilibrium stationary state, denoted (-)ncq, up to 
1/N 2 corrections. 

In particular, one may compute the average current (J2|) 
with respect to the measure ^ and get (with fa — T^ 1 ), 



(ji 



-viTJiTi+i-T^ + Oil/N 2 ), 



(7) 



where, i/(T<) = (5a(<7; - qi+i)\vl - pj+i\ + )T z is readily 
interpreted as the average collision frequency between the 
neighbors i and with respect to a global equilibrium 



measure at temperature Tj, Eq. ([5]). This computation 
therefore shows that the conductivity k(T,-), defined as 
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Ti 



(8) 



is identical to v{Tj). Furthermore a simple scaling argu- 
ment shows that k(T,) = v{Ti) = y/Tiv, where v denotes 
the collision frequency computed at unit temperature. 
Being the result of an equilibrium integration, the fre- 
quency may be computed with arbitrary precision. 

In order to get a better picture of the process that is 
described by the Boltzmann equation ([3j, we linearize 
the equation around the global equilibrium solution 
Doing so, we obtain an equation similar to ([3]), but with 
the collision operators L c ii±1 now replaced by L\ ln , 



j t fi(p,q,t) = -p- Vq/i + L^ 



2Lf»f, (9) 



where the linearized collision operator L' ln is obtained 
from Eq. (j4|) by replacing fa±i by equilibrium distribu- 
tions at common inverse temperature (3. 

The interpretation of the stochastic process described 
by the linearized collision operator is straightforward: 
when collisions take place, the particles velocities are up- 
dated as though they collided with stochastic thermal 
walls at inverse temperature (3 At each collision, the 
new velocities are independent from the previous ones. 

With this prescription, we now compute the conduc- 
tivity using the Green-Kubo formula, which is derived as 
follows. Integrated over time, the energy current between 
sites i and i + 1 takes the form 



1 vtMf] , 



Ji([0,t])= I ' t j i (s)da = i [P^ s ") 

J ° 0<s k <t 

(10) 

where the (sk) fceN are the successive collision times be- 
tween particles i and i + 1. The Green-Kubo formula, 
reads in our case, 



kgk(T) = 



1 1 N 
x lim — > 

2NT 2 t^oo t ^ 

i,k=l 



Ji([0,t])J k ([0,t]) 



(11) 

Using the expression (JTUJ) , @ , translation-invariance and 
the independence of the transfer of energy at each colli- 
sion, we get, after some calculations, 



1 

8T2 



(12) 



(Mqo - q0bo x -pi l + [(^) 2 - (pi) 2 } 2 )^ 



which, after further computations turns out to be equal 
to the collision frequency, kgk(T) — v(T). 

The square-strings model displayed in Fig. [T] lends 
itself to a detailed study of the dependence of the ratio 
njv on the parameter values a. 
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To this end we consider systems of varying sizes N 
with both ends in contact with stochastic thermal baths 
at respective temperatures XI = 1/2 and T + = 3/2. This 
gives rise to non-equilibrium stationary states with tem- 
perature profiles such as displayed in Fig. ^ which, as 
N increases, approach the corresponding solution of the 
heat equation, d x [n(T(x))d x T(x)] — 0, with k(T(x)) oc 
y/T(x). The ratio n/ v is obtained by linearly extrapo- 
lating to N — > oo finite N measurements of the spatial 
averages of n{Ti)/v{Ti), with n(T t ) defined by Eq. © 
and v{Ti) the collision frequency at the local tempera- 
ture, as functions of 1/N. 
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FIG. 2: (Color online) Non-equilibrium temperature profiles 
of the square-strings model with a = 2.08, for increasing val- 
ues of N — 5, 10, ... , 50. The black curve is the stationary 
solution of the heat equation. The inset displays the cor- 
responding measurements of (ft(Ti) / 'i/(T<)) . The infinite N 
extrapolation, k/v — 1.0037, is the approximate heat con- 
ductivity reported in Table [T] 

These values are reported in Table |TJ Notice the ex- 
cellent agreement with the prediction k = v, Eq. (fl~2f . 
as the value of the parameter a gets closer to its max- 
imal allowed value, the limit of rare collisions, in close 
agreement with the results presented in 0] for a class of 
coupled semi-dispersing billiards. In particular we un- 
derline that the parameter range of validity of our result 
is very similar to that observed in [j| , which further vali- 
dates that it is independent of the detailed nature of the 
local dynamics. 



%/5 — a 
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\/5 — a 
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a/5 — a 


k/v 


1.029 


6.000 


0.618 


1.892 


0.322 


1.320 


0.870 


3.159 


0.511 


1.632 


0.236 


1.0718 


0.736 


2.336 


0.413 


1.452 


0.155 


1.0037 



TABLE I: Measurements of k/v for selected values of a, ob- 
tained from data similar to Fig. [2] Our results indicate that 
c/i/-tlasa-> v 7 ^, in agreement with Eq. (JT2J ■ 

To summarize, we have showed that the derivation of 



Fourier's law in a large class of locally confined particle 
systems with hard-core interactions can be achieved from 
a Boltzmann-type approach with the main result that, in 
the appropriate limits, the heat conductivity is identified 
with the collision frequency. 

The same identity was derived in [8| in the context 
of semi-dispersing billiards. The comparison is interest- 
ing since, in contrast, chaos in the square-strings model 
results from a defocusing mechanism which takes place 
after particles interact. The identity between conductiv- 
ity and collision frequency therefore proves to be more 
general as it accounts for the transport properties of sys- 
tems lacking the local mixing property. In fact, the only 
dynamical property which is a priori necessary in our 
derivation is ergodicity of the local dynamics, i.e. in the 
absence of interactions. This property guarantees that 
two neighbors always interact provided the coupling is 
switched on, and that the fraction of time during which 
they interact is proportional to a fixed geometrical factor 
which can be adjusted by tuning the systems' parameters. 

The square-strings model is a perfect example of a sys- 
tem which lends itself with ease to a precise and reliable 
numerical analysis, while retaining the molecular chaos 
property. The square-strings model is actually a kind of 
higher dimensional fully chaotic stadium and displays a 
very rich structure of dynamical properties. 

We regard the proof of the molecular chaos hypothesis 
upon which our computation relies as a promising and re- 
alistic way to eventually obtain a clear picture of the dif- 
ferent mechanisms responsible for the origin of Fourier's 
law in a large class of mechanical systems. 
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